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Using 355 pb _1 of data collected by the CDF II detector in pp collisions at y/s = 1.96 TeV 
at the Fermilab Tevatron, we study the fully reconstructed hadronic decays B^ — > D^ir + and 
_B^ s j — > Dj" j7r + 7r + 7r~ . We present the first measurement of the ratio of branching fractions 

B{B° -» D~tt + tv + tt-)/B(B° -> D-tv + tt + tv~) = 1.05 ± 0.10(stat.) ± 0.22(syst.)- We also up- 
date our measurement of B{B° -» D~-k + )/B(B° -> D~ty + ) to 1.13 ± 0.08(stat.) ± 0.23(syst.) 
improving the statistical uncertainty by more than a factor of two. We find B(B® — ► D~ir + ) = 
[3.8±0.3(stat.)±1.3(syst.)] x 10" 3 and B(B° -> D^tt+tt+tO = [8.4±0.8(stat.) ±3.2(syst.)] x 1(T 3 . 

PACS numbers: 13.25.Hw 14.40Nd 



Hadronic B meson decays provide important informa- 
tion on both weak and hadronic interactions of heavy 
flavored mesons. The dominant hadronic decay modes 
of the B meson involve tree-level diagrams where the 
b — > c transition leads to a charmed meson and a vir- 



tual W boson, which often emerges as a charged ir, p, 
or ai(1260) meson The measurement of the ra- 

tios of branching fractions B(B° — > D~ir + {ir + n~])/ 
B(B° — * D~ir + [7r + 7r~]) reveals information about B 
decay mechanisms. One can attempt to separate the con- 
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tributions of various processes in B° — ► D~ir + decay 
and then predict B® — > D~ir + branching fraction using 
SU(3) 13,13) an d further estimate flavor SUY3) symmetry 
breaking effects which can be sizable @- 

Experimentally the sample of B® — > D;T7r + 7r + 7r~ 
decays is included in the data set used for the first 
measurement of the B® oscillation frequency Q. In 
this letter we present the first measurement of the 
ratio of branching fractions R(D3tt) = B(B® — > 
D~n + n + ir~)/B(B° — > Z?~7r + 7r + 7r~) using Dj" decays to 
</>7r~, K*°K~, 7r~7r~7r + , and decays to X + 7r~7r~. 
We also update our previous |8| measurement of the ra- 
tio R(Dn) = B(B° -> D~ir+)/B(B° -> £>-■*+). We 
measure the ratios of branching fractions because most 
of the systematic uncertainties cancel due to the simi- 
larity of final state kinematics. In addition, some com- 
mon theoretical factors [j| also cancel in taking the ratio. 
The measurement is performed using a sample of inclu- 
sive heavy flavor decays selected by the displaced track 
trigger described below, corresponding to an integrated 
luminosity of 355 pb _1 of pp collisions at y/s = 1.96 TeV. 

The components of the CDF II detector relevant for 
this analysis are briefly described below. A more com- 
plete description can be found elsewhere 0. Charged 
particle tracks are reconstructed in the pseudo-rapidity 
range \t]\ < 1.0, where rj is defined as — lntan(0/2), and 6 
represents the angle between the particle and the proton 
beam direction jll| . Tracks are reconstructed from hits 
in the silicon microstrip detector (SVX) and the central 
outer tracker (COT). Both detectors are inside a 1.4 T 
solenoidal magnetic field. The SVX detector is composed 
of LOO (single layer of silicon placed close to the beam 
pipe) , SVX II (five cylindrical layers of double-sided sen- 
sors), and ISL (outermost layer of silicon), providing up 
to 8 coordinate measurements in the r-cj) view |l2|. Sur- 
rounding the SVX is the COT, an open cell drift chamber 
with 96 layers of sense wires [13]. A sample rich in charm 
and beauty hadrons is selected by a three-level displaced 
track trigger. At level 1, tracks are reconstructed in the 
COT by the track trigger processor (XFT) 14]. The trig- 
ger requires two tracks with transverse momenta pt > 2 
GeV/c and the scalar sum pti + Pt?. > 4.0 GeV/c. The 
level 2 silicon vertex tracker (SVT) |l5j associates SVX II 
r-<p position measurements with XFT tracks, providing a 
precise measurement of the track impact parameter (do), 
i.e. the distance of closest approach of the track helix to 
the beam axis in the transverse plane. Decays of heavy 
flavor particles are identified by requiring two tracks with 
0.12 mm < do < 1 mm and an opening angle in the trans- 
verse plane 2° < \A<f>\ < 90°. A requirement L xy > 0.2 
mm is also applied, where L xy is defined as the distance in 
the transverse plane from the beam line to the two-track 
vertex projected onto the two-track momentum vector. 
The level 3 trigger receives the complete detector infor- 
mation and performs a full event reconstruction applying 
selection similar to level 1 and level 2 on offline quality 
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FIG. 1: Comparison of the sideband subtracted mass spec- 
trum of tt + tt + iy~ from B° and B® decays. Only the D~ — > 4>ty 
channel is used for B®. The B° histogram is normalized to 
B° 



quantities. 

B candidate reconstruction starts with a collection of 
tracks. No particle identification is explicitly used, and 
tracks are assumed to be either a pion or kaon to match 
the reconstruction hypothesis. A set of unique track com- 
binations making the <fi, K*°, D~ , Dj, and B candidates 
is formed. The track combinations reconstructed in 3 
dimensions must be consistent with forming a vertex, 
and combinations that fall outside a wide mass window 
around the mass of the respective meson are rejected. 

The Monte Carlo simulation is an essential part of this 
analysis. It is used to optimize the selection cuts, model 
signal and background, and to study the trigger and re- 
construction efficiency. We generate single B hadrons 
with the program BGenerator ITbI. The B-hadron de- 
cays are simulated with EvtGen [Tj| . This package has 
been extensively tuned by experiments at the T(45) res- 
onance and reflects the measured properties of B and D 
meson decays. The simulation is validated by comparing 
kinematic quantities of the decays and quantities relevant 
for the trigger selection with data. 

The selection requirements used to reject combinato- 
rial background are optimized by maximizing S/yS + B 
for each mode individually. The number of signal events 
(S) is derived from a Monte Carlo simulation of the 
CDF II detector and trigger. The number of background 
events (B) is estimated using data in the high-mass side- 
band interval, m(B) + 10 <j(B) to m(B) + 16 cr(B), where 
m(B) is the fitted mass and a(B) w 15 MeV/c 2 is the 
width of the signal peak. This sideband represents the 
combinatoric background underneath the signal peak. 
Selection requirements include cuts on the impact pa- 
rameter of the B meson, the Xr-d> °f the B vertex 



■5 



t rr 1 — e 

B s candidate Mass [GeV/c ] 



■ 'combinatorics 

L _■ 

[]B, -> D\ X 




combinatorics 



600 - 
400 - 
200 - 



j | B D -> D' n, D*p 






(f) 


i : : combinatorics 














| | B° -> D X 








^B'^D 2jiK, A b , B 


,-» D s 3n 




1st, d 


i ■ other physics 






"""" nT 







B s candidate Mass [GeV/c ] 



5 5 5 2 

B d candidate Mass [GeV/c ] 



FIG. 2: Mass spectra for (a) B° -> DJ% + (b) B° -> D; 
Dj{ir-TT-iv + )-K + TT + ir- (e) B° -» I)-7r + (f) 5° -> D-tt+TT+tT 
B -> DJX and B ^ D~ X decays. 



/)7r-)7r + 7r + 7r- (c) B° -> L»7(A"*°^-)7r + 7r+7r- (d) B° -> 
The "other physics" category corresponds to the inclusive 



fit in the transverse plane, the pt of the pion from the 
B decay in B® s ^ — > _D^7r + , and a minimum pt require- 
ment of the tracks for decays with 6 tracks in the final 
state. We exploit the narrow — > K + K~ resonance and 
A'* — > A"+7r~ resonance to suppress background by re- 
quiring 1010 MeV/c 2 < m(0) < 1029 MeV/c 2 and 840 
MeV/c 2 < m(K*°) < 940 MeV/c 2 . There are also re- 
quirements on L xy /a(L xy ) - the significance of the mea- 
surement of L xy for B and D vertices. 

The assumptions on the relative contributions of res- 
onant ai,p7r, and non-resonant 7r + 7r + 7r~ in the B° — > 



D~ 



signal affect R(D3tt) because Monte Carlo 



simulation shows that their reconstruction efficiencies dif- 
fer. We find that the contribution of the <Zi resonance is 
dominant. The 7r + 7r + 7r _ mass distributions were com- 
pared in data between B° and B® mesons and are com- 
patible within statistics, as shown in Fig^ The resonant 
fractions in B® — > D~ir + ir + Tr~ and B° — > D~-k + tt + it~ 
decays are assumed to be identical. 

To extract R(D3tt) (or equivalently R(DttY), we follow 
the general procedure described in Ref. [gj, using the 
following formula: 



R{D3ir) 



h <B°) B(D-) N(B° S ) 
/. ' c(Bj)' B(D7)' N(B°)' 



(1) 



where N(B®) and N(B°) are the measured signal yields, 
e(B )/e(B°) is the ratio of trigger and reconstruction ef- 
ficiencies extracted from Monte Carlo simulation, fd/f s 
is the ratio of b quark fragmentation fractions into B° 
and B° s mesons, and B(D~)/B(D~) is the ratio of the 



world average values for branching fractions of D and 
D~ mesons into the reconstructed final states [lij . 

The yields, N(B®) and N(B°), are extracted from 
the mass spectra in Fig|21 using a binned likelihood fit 
and are summarized in Table [I] In the normalization 
modes B° — ► D~ir + and B° — > D~ir + ir + ir~ we observe 
8098 ± 114(stat.) and 3288 ± 76(stat.) signal candidates 
respectively. The signal peaks are modeled with a sum of 
two Gaussians with the same mean values but different 
widths. We follow the treatment of backgrounds as done 
m Ref. 0. The combinatorial background is modeled 
with a sum of an exponential function and a constant. 
The shapes of other physics backgrounds are modeled us- 
ing Monte Carlo simulation, and their parameterization 
is fixed in the fits to the mass spectra. 

For B° meson candidates, the structure below the main 
B peak has two pronounced peaks. This structure is the 
result of a B° — > D*~7r + (7r + 7r + 7r~) decay followed by a 
D*~ — > D~ir° decay. Angular momentum conservation 
in B decay implies a spin alignment of the D*~ that, 
combined with the low Q value of of its decay, leads to 
the two distinct peaks in the £>~7r + (7r + 7r + 7r~) mass dis- 
tribution seen in the data. In the case of the analogous 
B° — > £>*~7T + 7r + 7r~ decay, the D*~ decays to Dj and 
a photon 94% of the time 0]. This produces a single, 
relatively broad structure as seen in our data. None of 
these backgrounds contribute under the signal peak. 

There are several backgrounds whose mass distribu- 
tions peak near the signal region and must be subtracted. 
They are the Cabibbo suppressed decays B^ — > D^K + 



6 



Decay Yield e(B°)/e(B°) B{p-)/B{Dj) /.//„ ■ BR{B°) / BR(B°) 

D-(<t>TY-)n + 494 ± 28 0.913 ± 0.004 4.40 ± 0.59 0.292 ± 0.020(stat) ± 0.0l2(syst) 

Dj ((f>TT~)TT + iT + iT- 160 ±17 0.814 ±0.010 4.40 ± 0.59 0.263 ± 0.029(stat) ± 0.018(syst) 

D- (K*K-)ty+tt + ty- 90 ± 17 0.352 ±0.009 3.80 ±0.77 0.274 ± 0.053(stat) ± 0.030(st/st) 

L> s "(7r-7r + 7r-)7r + 7r + 7r- 49 ± 11 0.397 ± 0.009 7.80 ± 1.50 0.293 ± 0.067(s£at) ± 0.021(st/st) 



TABLE I: Summary of event yields, ratios of efficiencies, and individual branching ratio measurements. The uncertainties 
listed on the yield and the ratio of efficiencies are statistical only. B(D~) /B(D~) is the ratio of B(D~ — » K + tt~-k~) to 
the corresponding branching fraction of the DJ" meson {B{D~ -> <tnr-,<j> -> K~K + ), B(D~ -» K*°K-,K*° -> K~ ~ir + ) or 
B(D~ -> -K-n+n-)). 



and B® s j — > D,^K + ~k 7r + , and mis-reconstructed baryon 
decays TVg — * A+7r~ and A° — > A+7r~7r + 7r~. The ra- 
tio of Cabibbo-suppressed B^ — ► D^K + background 
to the corresponding signal was fixed to the world av- 
erage ratio of branching fractions The ratio of 
Cabibbo-suppressed B®^ — > D^K + tt^tt + decay to the 
\v I 2 

signal is fixed to y"° 2 ~ 0.05. The fraction of Ab back- 
ground was fixed using the recent CDF measurement of 
B(A° b -> A+tt-)/6(B° -> Z3-7T+) [13. In all the cases, 
the ratios are corrected for the relative trigger and recon- 
struction efficiencies. 

In the case of the B® — > D~ir + ir + ir~ decay, there 
is a reflection from B° — > D~ir + ir + ir~ . If one of the 
pions from a D~ — > K + tt~it~ decay is reconstructed 
as a kaon, a peak is produced under the B® signal re- 
gion. These events contribute in the B® — > D~n + n + iT~ , 
Dj — > K*°K~ decay because the i4T*° resonance is 
broad. The fraction of B° — * D^tt + tt + tt~ events under 
B® — > _D77r + 7r + 7r _ peaks is calculated from the observed 
number of i?° mesons. The number of £?° — > Z)~7r + 7r + 7r~ 
events under the -Bl? — * Dj" (/^*°i ; ir _ )7r + 7r + 7r _ mass dis- 
tribution is estimated to be 141 ±6(stat.). The system- 
atic uncertainty assigned due to the B° — ► D~Tr + Tr + ir~ 
background subtraction is dominant in this channel and 
is a part of the fit model" (see Table |HJ|. For 
B® — > D,T7r + 7r + 7r~ with Z)~ — > 07r _ the corresponding 
B° background fraction is very small due to the narrow 
width of the (p. The contamination of the double charm 
B° — > D~D+, with Df — > Tr + ir~ir + , is estimated to 
be w 1% of the B — > Z?~7r + 7r + 7r _ signal and is sub- 
tracted from the measured yield. The contamination of 
B® — ► D~Df in the J3° signal is found to be negligi- 
ble. Applying a cut on the mass of three pions in the 
B® — > Z?77r + 7r + 7r _ decay around the mass of the Dj 
meson does not change the measured yields. 

Systematic uncertainties due to the fitting procedure 
are estimated by comparing the fitted yields after chang- 
ing the mass range in which the fit is performed and 
also by varying the parameters of functions describ- 
ing the backgrounds. The systematic uncertainty on 
the ratio of efficiencies e(B^)/e(B°) comes from various 
physics sources. In all cases, the systematics were esti- 
mated by observing the change in the ratio of efficien- 



cies when the effect was considered. The effect of the 
choice of B meson spectrum used by Monte Carlo sim- 
ulation was determined by reweighting the Monte Carlo 
events with px spectrum based on NLO calculations 0] 
to match the pt spectrum measured at CDF [2l|. To 
estimate the systematic uncertainty due to B and D life- 
times we varied the assumed lifetime of B and D me- 
son in signal Monte Carlo within world average values 
[l9j . The composition systematic applies to the decay 
B® — > D~ (ir~ir + ir~)ir + ir + iT~ only and is due to the 
limited knowledge of the resonances in D~ — > tt~tt + tt~ 
decay. We assign a systematic uncertainty due to the 
unknown resonance structure of the tt + it + tt~ system in 
B? . — > D^ir + ir + ir~ decay by varying the fraction of a\ 
component in both B° and B° s signal Monte Carlo in a 
range consistent with the observed shape. The system- 
atic uncertainties are summarized in Table [H] 

The results of the measurements for R(D3ir) are sum- 
marized in Table [I] To average the results of the mea- 
surements by the expected yield we use Eq.Q, where 
N(B®) is a sum of the yields in three B® channels and 
e(-Bg) is a linear combination of Monte Carlo efficien- 
cies multiplied by the ratio of the branching fraction of 
D~ decay in a given channel and B(D~ — > <pn). Using 
f s /f d = 0.259 ± 0.038 0, we obtain: 

R(D3n) = 1.05±0.10(stat.)±0.07(syst.)±0.14(br)±0.15(pr) 

The (br) and (pr) uncertainties refer to the uncertainty 
on the D meson branching fractions and the ratio of frag- 
mentation fractions f s /fd- 

Using Eq.Q and the input from Table HI we obtain: 

R(Dir) = 1.13±0.08(stat.)±0.05(syst.)±0.15(br)±0.17(pr) 

The measurement of R(Dn) is consistent with the pre- 
vious CDF measurement |£| and supersedes that result 
with a statistical uncertainty reduced by a factor of two. 
This measurement is consistent within uncertainties with 
the theoretical prediction of 1.05 ± 0.24 Q. 

In conclusion, we have presented the first measurement 
of the ratio of branching fractions R(D3ir). We also 
have measured the ratio of branching fractions R(Dir) 
improving the statistical uncertainty by more than a fac- 
tor of two. Using the world average values for B(B Q — > 
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D-TT+) and B{B° -> D-tt+tt+tt-) [19| we find 

) = [3.8±0.3(stat.)±1.3(syst.)*]xl0- 3 andi3(£° 
7T+7T-) = [8.4 ± 0.8(stat.) ± 3.2(syst.)] x 1CT 3 . 



7T 



Effect 


Syst. 
B -> Dtt 


Uncertainty [%] 

B -> D3iv 


B pr spectrum 


±3.0 


±3.0 


B° lifetime 


±2.1 


±2.1 


37r resonance structure 


n/a 


±2.5 


_D7 — > 37T composition 


n/a 


±3.0 


trigger simulation 


±1.2 


±1.1 


B° fit model 


±0.5 


±1.3 


B° fit model 


±1.5 


±(3.6-9.3) 


Total 


±4.2 


±(6.7-10.9) 



TABLE II: Summary of the relative systematics uncertainties 
on R(Dn) and R(D3tt). The range of values appearing in the 
second column reflects the differences in the fit systematic of 
the three B° decays. 
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